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c~| ' Masses and widths of the four hght scalar mesons a, k, ao(980) and /o(980) may be 

reproduced in a model where mesons scatter via a qq loop. A transition potential is used 
Q^l to couple mesons to qq at a radius of ~ 0.57 fm. Inside this radius, there is an infinite 

bare spectrum of confined qq states, for which a harmonic oscillator is chosen here. The 
coupled-channel system approximately reproduces the features of both light and heavy meson 
spectroscopy. The generation of a, k, ao(980) and /o(980) is a balance between attraction 
^ ■ due to the qq loop and suppression of the amplitudes at the Adler zeros. Phase shifts increase 

. more rapidly as the coupling constant to the mesons increases. This leads to resonant widths 

which decrease with increasing coupling constant — a characteristically non-perturbative 
effect. 

PACS numbers: 14.40. Cs, 14.40.Ev, 13.25.-k, 13.75.Lb 

Recent experiments have improved parameters of a, k, ao(980) and /o(980) [l]-[8], and there 
are now extensive data on their couphngs to decay channels. Their masses do not conform to the 
pattern set by well-known qq nonets, e.g. p, uj, K*{890) and 0. The mass of the ao(980) is ~ 450 
MeV above the mass of the a, unlike the near degeneracy between p and u. This has led Jaffe 
to suggest that they are 4-quark states [9]. However, Scadron has shown that the existence of a 
light scalar nonet may be a direct consequence of dynamical chiral-symmetry breaking [10]. 
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Here, we present a simple model which reproduces their essential features. It fits tttt and Ktt 
elastic phase shifts, and also the hne-shapes of ao(980) and /o(980), i.e., Uq — > nrj, KK and nrj', 
and /o — > TTTT, KK and rjr). 

The model in its two versions has been described in earlier publications [11] [12]. In one formu- 
lation [11], an explicit harmonic oscillator potential, with corresponding wave functions, is used 
for the bare qq states; it reproduces the gross features of the spectra of both light and heavy 
mesons. In the other approach [12], use is made of a so-called Resonance Spectrum Expansion, 
which allows, in principle, the use of any confinement spectrum for the bare states. In the present 
paper, the latter method is applied, though again with a harmonic oscillator. This choice is not 
crucial, as e.g. the /o(1370) lies well above states considered here and couples mostly to 47r. So 
the more conventional funnel potential is expected to give similar results. 

The confining potential is joined at radius tq to plane waves for meson pairs. Decays of qq 
states trapped in the confining potential are described by Fig. 1. The cr, /o(980) and ao(980) 
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Figure 1: Schematic picture of the model. 

appear as extra states through their dynamical coupling to this qq loop. If it were not for this 
coupling, they would be plane waves in the continuum. Coupling of mesons to the loops induces 
some degree of mixing of regular qq states with a and its relatives. 

A universal coupling constant A is used for quarks of all flavours, after a mass scaling as a 
function of the reduced quark mass ensuring approximate flavour blindness of our equations, 
provided the same scaling is applied to the radius Tq [13]. The relative 3-meson couplings at 
each vertex are taken from Tables I and II of Anisovich et al. [14]. These are the usual S'f/(3)- 
symmetric couplings for OZI-allowcd decays, to be contrasted with the totally flavour-symmetric 
choice used in the second paper of Ref. [11]. The latter convention gives rise to equally good flts, 
but results in a significantly larger value of A for the coupled a"-/o(980), as compared to the k and 
ao(980), in the present restriction to pseudoscalar-pseudoscalar (PP) decay channels only. The 
formulae allow for the pseudoscalar mixing angle Qps in the couplings of r] and 77', 

77 = nricos(©o + ^ps) - sssin(©o + ^ps) , 
rj' = nnsin(©o + ©Ps) + sscos(©o + ©Ps) , 

where ©q = arctan\/2 ^ 54.7°, thus defining Qps in the octet-singlet basis. The intrinsic (or 
bare) mixing between a and /o(980) is likewise expressed, this time in the nn-ss basis, in terms 
of a scalar mixing angle ©5 by 

(7 = nncos©^ — sssin©5 , 
/o = nnsin ©5 -|- sscos ©5 . 

Coupled relativistic Schrodinger equations are solved outside and inside the radius Tq, and are 
matched at that radius. The resulting multichannel T matrix can be written down in closed form. 
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In the case of the k and the ao(980), as well as the decoupled a and /o(980), there is only one qq 
channel, resulting in the relatively simple expression 



(3) 



n=l Ks=0 



E-E. 



where A is the overall coupling, Tq is the delta-shell radius, gi{s) is the relative coupling of the ith 
(z = l,Nf ) two-meson channel, depending on the radial excitation (see Ref. [12], Table 1 of 4th 
paper). Eg are the levels of the discrete confinement spectrum, /Xj = fii{E) and ki = ki{E) are the 
relativistic reduced mass and momentum in the ith channel, and ji^ and h^.^ are spherical Bessel 
and Hankel functions, respectively. For the coupled a and /o(980), there are two qq channels, 
i.e., nn and ss, giving rise to a more comphcated expression. Writing Tij{E) — Nij{E)/D{E), we 
have 
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where the index a = 1, 2 refers to the nn and ss confinement channels, respectively. 

The formulae of Anisovich et al. allow for separate coupling constants for ss and nn; we find 
this refinement does not give any significant improvement, so we use the same value of A for both. 

A key feature of the model lies in the automatic inclusion of Adler zeros in scattering ampli- 
tudes, which results from using relativistic reduced masses for two-meson channels [15]. It leads 
to Adler zeros at s = for tttt and KK, at s = — nri^ for rj-n, and at s = m\ — for Kn] 
although these are not quite conventional values, e.g. 0.5 for tttt, the difference is insignificant 
for present purposes. The Adler zeros are essential to reproduce tttt and Xtt elastic phase shifts 
near threshold [16]. 

Experimental data are used to determine the coupling constant A and transition radius rg. 
The optimum values turn out to be in rough agreement with those in previous work [12] [13] (see 
below). Phase shifts for tttt nn are taken from the BES publication [2] on the a pole (Method 
1), where a combined fit was made to (i) BES data on J/^^ — > a;7r+7r~; (ii) Cern- Munich data [17]; 
(iii) Ke4 data of Pislak et al. [18]; (iv) the nn scattering length oq = 0.220 ± 0.005 determined 
from chiral perturbation theory (ChPT) by Colangelo, Gasser and Leutwyler [19]. Phase shifts 
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for Kn — > Ktt are taken from the combined fit [7] to LASS data for Kn elastic scattering [20], and 
E791 [4] and BES data [5] [6] for the k pole; because of uncertainties in fitting Kq{1A3Q) (discussed 
below), the fit is made from threshold to 1.2 GeV after subtracting phases due to i^^^(1430). The 
hne-shape of ao(980) is taken from the recent combined analysis [21] of KLOE data for (p 700 
[22], and Crystal Barrel data on pp —>■ cuao and oott [23]. The line-shape of /o(980) is taken from 
BES data [8]. 

The coupling constant A and radius ro are fitted for each scalar resonance separately, except 
for the /o(980). We also make a combined fit to cr + /o(980), in which both states emerge 
dynamically, by coupling the bare isoscalar nn and ss channels to vrvr, KK, 7777, rji]' and ri'r]'. 
Note that, although there is no direct coupling of ss to tttt, the higher-order OZI-allowcd process 
ss —>■ KK — > nn — > TTTT does produce a coupling after all. This also produces a dynamical 
mixing of nn and ss, in addition to the intrinsic quark-level mixing assumed above. Results for 
parameters and pole positions are shown in Table 1. Striking is the considerable effect on the 



Resonance 


A 


To 


Pole Positions 




(GeV-3/2) 


(GeV-i) 


(MeV) 


a 


2.92 


2.84 


555 - 1262 


K 


2.97 


3.29 


745 - 1316 


ao(980) 


2.62 


2.81 


1021 - 147 


(T + /o(980) 


3.11 


2.71 


530 - 1226 , 1007 - 138 



Table 1: Fitted parameters. For ao and /o, second-sheet poles are quoted. 

(7 pole from coupling to the /o(980), showing the influence of nn-ss mixing via the kaon loop. 
As for the parameters, there is a reasonable agreement for all cases. The A parameter varies by 
±9% and the radius parameter by ±10%. Moreover, both parameters are roughly compatible 
with the values A = 0.75 GeV~^/^ and ro = 3.2 GeV~^ used in previous model calculations with 
only one continuum channel [12] [13], taking into account that in the present study the squared 
relative couphngs of the included PP channels are normahsed so as to add up to 1/16, whereas in 
the single-channel case this number was just set to 1. Some small variations of parameters are to 
be expected due to the choice of a delta shell for the transition potential, which is very sensitive 
to the precise radial wave functions of the different PP channels. We shall comment on global 
features first and return to detail later. 

In Fig. 2(a), the dashed curve shows our fit to the BES fit (full curve) to tttt phase shifts 
after subtracting the /o(980) component. Structure in our fit near the KK and r]r] thresholds 
originates from a — > KK and rjrj. These thresholds were not included in the analysis of BES 
data, so the comparison illustrates the magnitude of possible effects from those thresholds. There 
is also structure in our fit from 800 to 1000 MeV, caused by quantum-mechanical mixing between 
/o(980) and the a amplitude. Figure 2(b) shows the fit (dashed) to Kt^ phase shifts up to 1.2 
GeV from LASS (full curve) after subtracting the component due to Kq{1A?,{)). Note that our 
model predicts a small cusp effect at the Krj threshold (also sec Ref. [24]), but within the errors 
of the data (typically ±3°). Figures 2(c) and 2(d) show fits to the line-shapes of /o(980) and 
ao(980), respectively. The original fits [21] to experimental data were made using Flatte formulae 
with widths of the form g'^p for each channel, where p is phase space. Our model is then fitted 
to these parametrisations, giving slightly more complicated line-shapes. Also notice the cusp at 
the rj'-K threshold in Fig. 2(d) due to the opening of that channel. 

Resonance formation arises from rescattering processes within the central loop of Fig. 1, i.e.. 
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Figure 2: Fits to (a) the o component of the tttt phase shift for elastic scattering; (b) the k 
component of Ktx phase shifts for elastic scattering; (c) + (d) the line shapes of /o(980) and 
ao(980) in tttt resp. 7777 channels. 



unitarisation effects. In this respect, the model has some similarity with the approach of Oset, 
Oiler, Pelaez et al. [25]-[30], who take their meson-meson Born terms from ChPT, and then 
unitarise the amplitudes (see Ref. [31] for an alternative unitarisation scheme). Note, however, 
that the present model is unitary from the start, via a coupled-channel (relativised) Schrodinger 
equation. Moreover, explicit quark degrees of freedom with a complete confinement spectrum 
are present, contrary to ChPT. We also observe that Ref. [28] requires a "pre-existing" s-channel 
flavour-singlet resonance at 1021 McV in order to reproduce the /o(980). In our model, such 
a low-lying s-channel contribution is not needed, only a regular bare P-wave ss spectrum from 
~ 1.5 GeV upwards. The /o(980) then appears spontaneously, though with a tttt width which 
depends on the bare scalar mixing angle 65 as well as the dynamical mixing with a. The model 
should also be contrasted with meson-exchange models (see e.g. the work of the Jiilich group, 
Janssen et al. [32]), in which attraction is generated by a few t-channel meson exchanges, whereas 
in our approach it arises from the coupling to an infinity of bare s-channel qq states. Furthermore, 
light scalar mesons dynamically generated in the Jiilich model need not show up as a complete 
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nonet. Nevertheless, the two models share the feature that the stronger the couplings, albeit of 
different natures, the narrower the resonances (see also below). Most akin to our model is the 
one of Tornqvist [33], also based on bare s-channel states coupled to free mesons, which however 
does not predict the k meson. 

A key point is that phase shifts rise more rapidly as the coupling constant A increases. Con- 
sequently, the masses of the low-lying resonances move with A, but their widths decrease as A 
increases. This is illustrated for the 4 light scalars in Table 2. We see that for very large cou- 



A 




a 


K 


/o(980) 


ao^yoU j 


1.5 


942 


- 1794 










2.0 


798 


- 1507 










2.2 


738 


- 1429 


791 - 


1545 




1081 - 18.0 


2.4 


682 


- 1368 


778 - 


1472 




1051 - i25 


2.6 


633 


- 1319 


766 - 


1409 


1041 - 113 


1024 - 145 


2.8 


589 


- 1278 


754 - 


1355 


1028 - 126 


998 - 161 


3.0 


549 


- 1243 


743 - 


1309 


1015 - 135 


978 - 160 


3.5 


468 


- 1174 


717- 


1219 


976 - 137 


896 - il42 


4.0 


404 


- 1123 


693 - 


il55 


948 - i38 


802 - 1103 


5.0 


308 


- 150 


651 - 


169 


889 - 134 


711 - 140 


7.5 


216 


+ 10 


610 - 


f 10 


752 - 125 


632 + 10 


10.0 


142 + iO 


560 -1- 10 


633 - 117 


577 + 10 



Table 2: Movement of the cr, k, /o(980) and ao(980) poles as the coupling constant A Is varied. 
Bound states are Indicated by "-I-10" . Units are MeV for the poles and GeV~^/^ for A. 

pling, the scalars tend to become bound states (also see Ref. [12], 3rd paper), while for small 
coupling they either disappear to —ioo in the complex energy plane or become virtual states. 
In this respect, o", k, /o(980) and ao(980) behave in a completely different fashion to regular qq 
nonets, (e.g. p, i^*(980), etc.) whose resonance widths are directly proportional to coupling 
constants. One of us attempted to fit 10 ratios of coupling constants for cr, k, ao and /o to this 
conventional scheme [21]. This attempt failed to reproduce many observed ratios, particularly 
for 5'^(/t K-n)! g^io tttt) which came out a factor 3 smaller than experiment. The reason for 
the failure is clear: non-pcrturbative effects due to the qq loop are crucial. In the present model, 
the ratio of k and a widths is reproduced using closely similar A parameters. 

Studying in more detail the trajectory of e.g. the a pole, we see how it emerges from the 
scattering continuum, passes the optimised position in Tablel for A ~ 3 GeV~^/^, and then 
approaches the real axis. However, the pole clearly slows down for increasing coupling, and 
very large values of A are needed to make a a tttt bound state. Here we see the Adler zero in 
action, whose influence increases as the pole moves to lower energies. The k pole behaves in 
a similar way, though it is more difficult to trace for small A because of the Kvl Adler zero at 
~ 800 MeV. Descotcs-Genon and Moussallam find a k mass of 658 ± 13 MeV and a width of 
r = 557 ± 24 MeV [34]. Their low mass does not fit well with what we find in Table 2. We 
remark that their calculation docs not include dynamics of the strong coupling of Kr( to K-n and 
consequent mixing between ii'Q(1430) and k. Below the Kvl threshold, the analytic continuation 
of that channel definitely plays a role, as our calculation verifies. The Kvl Adler zero could affect 
details of their calculation, and their errors may therefore be too small. 

The role of the Adler zero is crucial in understanding cr, k and Cq. In particular, the Adler 
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zero for oq prevents it from having a mass close to the tttj threshold. This is why it settles close 
to the KK threshold, since that channel has only a distant Adler zero. We have examined the 
pole structure of the nr] amplitude. It is interesting that, as well as having a narrow second sheet 
pole due to ao(980), there is a very broad pole at 991 — i815 MeV, resembling a bit the broad 
poles due to a and k; this broad structure is too wide to determine experimentally. 

It is important to test the stability of the fits to data. Of many such tests, one is to examine 
the stabihty of pole positions over the complete range of parameters shown in Table 1. This leads 
to the range of pole positions shown in Table 3. The poles do not disappear for any parameters 



Resonance 


Range of Pole 




Positions (MeV) 


(7 


(476-628) - 1(226-346) 


K 


(738-902) - 1(286-434) 


/o(980) 


(989-1040) - i(ll-38) 


ao(980) 


(960-1049) - 1(23-83) 



Table 3: Range of pole positions for the complete range of A and ro parameters of Table 1; 
second-sheet poles are quoted for /o (only as coupled to a) and Gq. 

within this range. We have already seen that the a and k poles even survive for a much wider 
range of A values alone. So their existence is secure. However, it is a different story for /o(980). 
If A is reduced to 2.7 GeV-^/'^. the cusp at the KK threshold begins to soften, and at Ri2.5 the 
pole breaks away from the KK threshold to become a virtual state at ~ 1050 MeV. 

The pole position and line-shape of ao are likewise sensitive to A. As it is decreased, the cusp 
at the KK threshold begins to soften for A = 2.3 GeV~^/^. The ao pole breaks away from the 
KK threshold and becomes a virtual state at about 1090 MeV for A ~ 2.1. It is clear that the /o 
and ao are less tightly bound than a and k. Jansscn et al. [32] have noted that ao(980) appears 
as a virtual state in their calculation based on meson exchanges. 

The model is undoubtedly an over-simplification; this is dehberate, so as to expose the essential 
features. However, these simplifications lead to difficulties with some details, which we now 
discuss. Firstly, one of the major unknowns afflicting all spectroscopy is how to continue effects 
of inelastic channels below their thresholds, e.g. in the processes KK tttt and Kr]' — > Ktt. 
This is relevant to details of both Kir and tttt phase shifts, particularly the former. The model 
predicts a large coupling of Kq(1A30) to Krj', in agreement with LASS data [20]. However, it 
also predicts a similar large coupling of k to Krj'; this is because mesons couple via the qq loops 
of Fig. 1. The combined effect is to predict an even larger amplitude for Kr]' — > Kir than is 
fitted experimentally. It is important to attentuatc this amplitude below the Kr]' threshold, to 
avoid distorting the fit to the k. The model builds the Ktt Adler zero into the amplitude for 
Krj' Ktt. However, the fit improves if we introduce an additional form factor exp(— 4|A;|^), 
where k is the magnitude of the kaon momentum below the inelastic threshold in GeV/c. This 
corresponds to the fact that ^0(1430) has a long tail of Kr]' near its threshold. For consistency, 
the same suppression factor of subthreshold contributions is applied to the other scalars, too. 

A further technicality is that, coupling only PP channels, too narrow a width is predicted for 
K*{U30y. ~ 200 MeV, compared with 294 ± 23 MeV quoted by the PDG [35]. To avoid this 
conflict, we take the k amplitude from the flt to experimental data in Ref. [6], and flt it only up 
to 1.2 GeV, where the effect of Kq{1430) in our model becomes negligible. Preliminary results 
indicate that the inclusion of additional inelastic channels, such as vector- vector and scalar-scalar. 
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may improve the widths of higher resonances. This will be the subject of future work. 

Next, the assumption of a sharp radius for the transition potential is an approximation; this 
affects radial wave functions. Thirdly, mesons outside the transition radius are treated as plane 
waves, but in reality they will be affected by i-channel meson exchanges, which for simplicity are 
absent in our approach. 

We have fitted a alone and /o(980) alone, but also allowed both to couple simultaneously to 
TTTT, KK^ rjr], r]r]' and r]'r]'. There is obvious mixing between a and /o(980). When a and /o(980) 
are fitted together, the optimum fit gives a bare scalar mixing angle of ~4°. However, if /o(980) 
is fitted alone with this mixing angle, the nn width comes out far too small, ~2 MeV, compared 
with the observed full width at half maximum of 34 ± 4 MeV. Clearly, the additional dynamical 
mixing of /o(980) with a, mostly through intermediate KK, is crucial for a realistic description. 

With the present sharp transition radius, it is difficult to fit simultaneously: (i) the elasticity 
parameter r) above 1 GeV; (ii) the intensity of tttt —>■ KK from 1 to 1.2 GeV; (iii) tttt phase 
shifts from 1 to 1.2 GeV. Any may be fitted alone, but the threefold combination exposes some 
confiicts. Figure 3 shows two fits. In (a), we show the recent tttt phases of Kaminski, Pelaez and 




M (GeV) M (GeV) 

Figure 3: Fits to (a) the nn phase shift with ro, A and 9s fitted to these data alone; (b) the 
elasticity parameter rj for tttt elastic scattering; (c) the magnitudes squared of the T-matrix for 
TTTT — > KK; and (d) as (a) with tq, A and ©5 fitted to all data simultaneously. 

Yndurain [36] above 1 GeV. Figure 3(a) and full fines on (b) and (c) show fits to individual sets of 
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data, varying A, ro and 6s- This achieves a good fit to all sets of data. Figure 3(d) and the dashed 
lines on Figs. 3(b) and (c) show the best fit when all three sets of data are fitted simultaneously. 
The pole position of /o(980) is then at 1007 — i38 MeV, compared with the experimental value 
(998 ± 4) - i(18 ± 4) from BES parameters. 

The optimum value of ©5 is ~4°, quite a bit smaller than frequently quoted values (see e.g. 
Ref. [37]). However, fitting /o(980) alone gives rise to a much larger ©s, in excess of 30°. So there 
is a sizable, s-dependent dynamical contribution to the scalar mixing angle, which is largest in 
the vicinity of the /o(980), that is, near the KK threshold. This is a physically appeahng result 
of our fits in this sector. 

A further detail concerns the pscudoscalar mixing angle Qps- The experimental value from 
Crystal Barrel is (— 17.3± 1.8)° [38]; that is the value used here (also see Ref. [24]). However, the 
fit to the line-shape of ao(980) — > t^tt is sensitive to Qps- For those data, the optimum value is 
Qpg = -(7 ± 3)°. This value is used for Fig. 2(b). With Qps = -17.3°, the fit is 15% too high 
on both sides of the peak, as shown by the dotted curve. These discrepancies might be associated 
with the extensive KK cloud around ao(980) or with the qqqq correlations proposed by Jaffe. 
Alternatively, small yet non-negligible contributions from higher, closed thresholds may improve 
the predicted Qps, as indicated by preliminary fits. 

Despite the blemishes described above, the essential results are as follows. The model of Fig. 1 
is capable of generating all of a, k, ao(980) and /o(980) dynamically via the intermediate qq loop. 
This is a non-pcrturbative effect. The overall coupling constant required for the four resonances 
varies by only ±9% and the radius parameter tq by ±10%. Moreover, all scalar poles survive, 
with reasonable values, over the complete range of fitted parameters, so the model is very robust. 
The widths of these dynamical resonances decrease as the coupling constant increases. This is an 
essential difference from regular qq states, where, for not too large coupling, widths are roughly 
proportional to the square of the coupling. Nevertheless, the formation of these four resonances 
cannot been dissociated from the presence of a bare qq spectrum, which is crucial for the generation 
of dynamical resonances with moderate widths. In particular, for each ground-state bare scalar 
qq state, a pair of resonances is produced by coupling to the meson-meson continuum, namely one 
light, non-standard scalar meson, and one regular scalar in the energy region 1.3-1.5 GeV. This 
phenomenon was first observed two decades ago (see 2nd paper of Ref. [11]), and later confirmed, 
to some extent, by Tornqvist [33]. 
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